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ABSTRACT 

This paper presents the thkd extension to the Very Large Baseline Array (VLBA) Calibrator Survey, con- 
taining 360 new sources not previously observed with very long baseline interferometry (VLBI). The survey, 
based on three 24 hour VLBA observing sessions, fills the areas on the sky above declination -45° where the 
calibrator density is less than one source within a 4° radius disk at any given direction. The positions were 
derived from astrometric analysis of the group delays determined at 2.3 and 8.6 GHz frequency bands using the 
Calc/Solve software package. The VCS3 catalogue of source positions, plots of correlated flux density versus 
the length of projected baseline, contour plots and fits files of naturally weighted CLEAN images as well as 
calibrated visibility function files are available on the Web at http : / /gemini .gsfc.nasa. gov/vcs3 
Subject headings: astrometry — reference frames — astronomical data bases: catalogues — quasars: general 
— radio continuum: general — reference systems: techniques — interferometric 



1. INTRODUCTION 

\ Catalogues of positions of compact radio sources with the 
highest precision are the basis for many applications. Among 
of them are imaging faint radio sources in phase referenc- 
ing mode, monitoring the Earth's rotation, space geodesy, 
and space navigation. The method of VLBI first proposed 
by Matveenko et al. ( 1965) allows us to determine positions 
of sources with nanoradian precision (1 nrad « 0.2 mas). 
The first ca talogue of so urce positions determined with VLBI 

' (jCohen & Shaffedl 19711) contained only 35 sources. It is de- 
sirable for the above mentioned applications to have many 
more sources. In 1998 the ICRF catalogue of 667 sources 
produced by analyzing VLBI observations made in the frame- 
work of space geodesy programs was published (MaetaL 
[1598). Six years later, the extension of this c atalogue, ICRF - 
Ext2, with 109 more sources, was published ( Fev et al.l2004h . 
However, the sky coverage of these catalogues was still not 
sufficient for some applications. Successful phase referenc- 
ing requires a calibrator within 4° from the target source with 
a precise position and known source structure. The proba- 
bility of finding a calibrator from the ICRF-Ext2 catalogue 
within 4° of any random target above -45° declination is 
about 50%. Also, the ICRF catalogue does not provide source 
maps. In order to overcome these deficiencies, twelve 24 
hour observing sessions with the VLBA, called VLBA Cal- 
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ibrator Surveys, were performed. Analysis of these obser- 
yations produce d the VCSl (jBeaslev et al. 2002) and VCS2 
JFomalontetal.ll2()0l catalogues of 1611 new sources not 
listed in the ICRF catalogue. Approximately 80% of these 
sources are suitable as phase referencing calibrators and as 
target sources for space geodesy programs. 

Combining all suitable calibrators from the ICRF-Ext2, 
VCSl and VCS2 catalogues we get a list of 1986 sources. 
Over 89.5% of the sky north of -45° declination we find at 
least one calibrator within a 4° radius disk at any given direc- 
tion. In this paper we present an extension of the VCS 1 and 
VCS2 catalogues, called the VCS3 catalogue, mainly concen- 
trating on the other 10.5% of the sky where the source density 
is lowest. Since the observations, calibrations, astrometric so- 
lutions and imaging are similar to that of VCSl, most of the 
details are described by Beaslev et al. (2002). In section 2 we 
describe the strategy for selecting the 450 candidate sources 
observed in three 24 hour sessions with the VLBA. In section 
3 we briefly outline the observations and data analysis. We 
present the catalogue in section 4, and summarize our results 
in section 5. 

2. SOURCE SELECTION 

Having unlimited resources one could try to observe all 
sources stronger than some limiting flux. However, only 
sources with bright compact components can be detected with 
VLBI and might be useful fo r phase refer encing or geodetic 
applications. IKellermann. Paulinv-Toth. & DavisI ( 1968) first 
showed that the distribution of sources over spectral index 
a (F oa v"') has two peaks: one near a = -1 (steep spectrum) 
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and another near a = (flat spectrum). Later, it was con- 
firmed that extended objects dominate in the steep spect rum 
population (see the review bv lKeUermann & Owe iHQB S*). In 
compact regions, which have the synchrotron mechanism of 
emission, the peak in the spectrum caused by synchrotron 
self-absorption has frequencies higher than 1 GHz due to their 
small size. Thus, if a dominating emission comes from com- 
pact regions, the spectrum of the total flux density will be 
predominantly flat or inverted. 

Examining spectral indexes of the PKScat90 catalogue 
(fWright and Olnincek 1990 ), which is believed to be com- 
plete at the 95% level for sources with flux density 
F > 0.25 Jy at 2.7 GHz jRolton. Savage. & Wright' TgrP). 
we found t hat only 26% of so urces have spectral indexes 
a > - 0.5. Mversetal. (2003) used NVSS (C ondon et alJ 
Il998h and GB6 (Gregorv et al. 1996) catalogues for selection 
of sources with spectral indexes > -0.5 at 5 GHz, and they 
found 1 1685 sources brighter than 30 mJy, or 24% of the total 
number of sources. It should be noted that a significant frac- 
tion of the sources with spectral index in the range [-0.5, -0.4] 
belongs to the steep spectrum population and are expected to 
have extended structure. Thus, taking sources at random, the 
probability that a source will be compact, and therefore, de- 
tectable with VLBI does not exceed 20%. A carefully de- 
signed strategy for source selection can significantly increase 
the yield of detections. 

We considered all bright sources at X-band (8.6 GHz) 
and S-band (2.3 GHz) with spectral indexes greater than 
-0.5 as candidates for the VCS3 campaign. In order to 
find these sources, we first analyzed catalogues of com- 
pact flat-spectrum sources: JVAS (Wrobeletal. 1998), 
CLASS (Mvers et al. 2003), Parkes qu arter- Jansky catalogue 
djackson et aL,2002) . VLA South (.Winn. Patnaik. & Wrobell 
l2003li " and selected all sources brighter than 100 mJy at decli- 
nations above -45° with arc distances from any known cali- 
brator greater than 3.9°. This produced a list of 374 sources. 
Not all these catalogues provide flux de nsity at more than one 
band. We used the CATS database ( Verkhodanov et alJl9'97h 
for cross-checking the reported flux density and for acquir- 
ing information about flux measurements for these objects 
from other sources. Unfortunately, most of these measure- 
ments were done at different epochs. Since many flat spec- 
trum sources are variable, the estimates of the spectral in- 
dex made at different epochs may be unreliable. Instanta- 
neous 1-22 GHz broad-band spectra at 4-6 frequencies for 
3000 AGN with flux density greater than 100 mJy at 2.3 GHz 
and a > -0.5, including many sources from this list, were ob- 
served with the 600 meter ring radio telescope of the Rus- 
sian Academy of Sciences RATAN-600 in tra nsit mode (see 
method of observations and data reduction in iKovalev et al.l 
Il999h as a part of an ongoing survey program. All flat spec- 
trum sources from this sample with flux densities greater than 
200 mJy at X and/or S band and not previously observed with 
VLBI, a total of 1 16 sources, were included in the VCS3 can- 
didates list. After cross checking we removed sources from 
the list for which the estimate of the spectral index was ei- 
ther unreliable or less than -0.5. The total number of sources 
selected from these catalogues was 300. 

After picking up these 300 sources from all known com- 
pact flat-spectrum catalogues we still had areas where the 
density of calibrators was less than the targeted density. In 
order to fill these areas we used the NVSS catalogue of about 
1.8 million sources ( Condon et al . 1998). First, we selected 
2704 sources with total flux density greater than 100 mJy at 



1 .4 GHz in the areas where the distance to the closest calibra- 
tor and already selected candidates was greater than 4°. The 
NVSS catalogue does not provide information about spectral 
index, so we checked these 2704 sources against the CATS 
database, which nowadays contains one of the most complete 
collections of radio catalogues, and found the sources ob- 
served under other campaigns at different frequency bands. 
If they were not found in any other catalogue, for example, 
they were absent in the GB6 catalog ue with flux limit 18 mJy 
and in the PMN (Wright et al. 1996) catalogue with flux Hmit 
40-70 mJy at 5 GHz, this indicates that with a high probabil- 
ity their spectral index is significantly less than -0.5. In total, 
this initial list had 302 sources. 

We examined the spectrum provided by CATS and subjec- 
tively ascribed a class to each source: 

— the source has flat spectrum, and the estimate of the 

spectral index is reliable (e.g., RATAN-600 instanta- 
neous spectrum is available), and the estimated corre- 
lated flux density is greater than 100 mJy; 

1 — the estimate of the spectral index is less reliable, or the 

spectral index is in the range [-0.4, -0.2], and the esti- 
mated correlated flux density is greater than 100 mJy; 

2 — the estimate of the spectral index is unreliable or in the 

range of [-0.5, -0.4], and/or the estimated correlated 
flux density is less than 100 mJy. 

Then for each source we computed a score based on angular 
distance to the closest calibrator, flux density, spectral index 
and class: 

S = 10' {{D-Du,n)/D,i,„)' + 105(a + 0.5)^ + 

2 • 10^(2-0 + IgF (1) 

where D — angular distance to the closest calibrator or al- 
ready selected candidate, Dum minimal angular distance (4°), 
a — spectral index, C — class, F — flux density in Jy ex- 
trapolated to 8.6 GHz. An iterative procedure computed the 
score, selected the source with the largest score and ran again. 
Using this method we selected 150 sources with the highest 
scores in addition to the 300 sources previously chosen. 

3. OBSERVATIONS AND DATA ANALYSIS 

The VCS3 observations were carried out in three 24 hour 
observing sessions with the VLBA on 2004 April 30, 2004 
May 08, and 2004 May 27. Since flie correlated flux density 
for selected candidates was expected to be different, we split 
our sources into three categories: 

a — sources of class brighter than 200 mJy; 

b — sources of class with flux density in the range 100- 
200 mJy or sources of class 1 brighter than 200 mJy; 

C — all other sources. 

Integration time was chosen to be 90 seconds for sources 
from category a, 180 seconds from category b, and 300 sec- 
onds from category C. The sources above -20° declination 
were scheduled for two scans at different hour angles in order 
to improve the uv coverage, while the sources in [-45°, 20°] 
declination range were scheduled only for one scan. The tar- 
get sources were observed in a sequence designed to mini- 
mize loss of time from antenna slewing. In addition to the 
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target sources, 77 strong sources were taken from the GSFC 
astrometric and geodetic catalogue 2004a_astro ^. Observa- 
tions of 3^ strong sources from this list were made every 
1-1.5 hours during observing sessions. These observations 
were scheduled in such a way, that at each station at least one 
of these sources was observed at an elevation angle less than 
20°, and at least one at an elevation angle greater than 50°. 
The purpose of these observations was to provide calibration 
for mismodelled troposphere path delays and to tie the VCS3 
source positions to the ICRF catalogue. On average, the an- 
tennas were on-source 70% of the time. 

The observations used the VLBA dual-frequency geodetic 
mode, observing simultaneously at 2.3 GHz and 8.6 GHz. 
Each band was separated into four 8 MHz channels (IFs) 
which spanned 140 MHz at 2.3 GHz and 490 MHz at 8.6 GHz, 
in order to provide precise measurements of group delays for 
astrometric processing. Since the a priori coordinates of many 
candidates were expected to have errors of up to 30", the data 
were correlated with an accumulation period of 1 second in 
64 frequency channels in order to provide extra-wide win- 
dows for fringe searching. 

Processing of the VLBA correlator output to obtain the as- 
trometric observables was made using t he Ast ronomical Im- 
age Processing System (AJPS) (Greisenl ll988h . First, phases 
of the calibration signals injected into the receivers were sub- 
tracted from the fringe phases for each scan and each chan- 
nel. Then the phases of a single strong calibrator scan were 
subtracted from each scan to flatten out the phases across the 
band of each individual channel. Sampler bias correction to- 
gether with initial amplitude calibration was performed using 
the recorded system temperatures and gain tables (tasks 'AC- 
COR" and "APCAL"). Task "FRING" was then run for each 
baseline separately (with SNR cutoff = 3) to obtain the resid- 
ual narrow-band group delays for each channel, residual sin- 
gle channel fringe phases, and residual single channel phase 
delay rates. Task "MBDLY" was then used to compute the 
residual group delays from the 4 single channel fringe phases 
at each frequency band. Then task "CL2HF" was used to 
compute the total group delays, phase delay rates and fringe 
phases by combining the residuals with the correlator model, 
and converting them from geocenter to reference station quan- 
tities. Finally, the AlPS output was reformatted for input into 
the Calc/Solve geodetic/astrometric VLBl analysis software 
package-' using task "HF2SV". 

The software package Calc/Solve was used for determina- 
tion of accurate positions of all detected sources. Astromet- 
ric analysis was performed in several steps. First, we used 
narrow-band group delays determined over each 8 MHz wide 
channel and incoherently averaged separately over all chan- 
nels within X and S band. Narrow-band group delays are 
less precise than wide-band group delays, but have very wide 
ambiguity spacings: 4 mks. In the first step X-band single- 
band delay data (and S-band data for the sources detected 
only at S-band) were used for estimation of new source po- 
sitions and clock functions modeled by second degree poly- 
nomials. These position estimates were used in an S-band 
wide-band group delay solution for short baselines between 
stations Pietown, Kitt Peak, Owens Valley, Los Alamos and 
Fort Davis in order to resolve the group delay ambiguities 
with spacings of 100 ns and then again solve for source po- 
sitions in the second step. In step three these estimates were 
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used in an S-band group delay solution for longer baselines 
for further ambiguity resolution. In the fourth step the previ- 
ous estimates of source positions were used for ambiguity res- 
olution of X-band group delays (with spacings 29 ns) on short 
baselines, and a new solution was made. These estimates then 
were used for resolving X-band group delays at long base- 
lines in the next step. In the sixth step ionosphere free com- 
binations of X and S band wide-band group delays were used 
for data analysis. The estimation model at this step included 
coefficients of linear splines for clocks of all stations, except 
a reference station, and troposphere zenith path delays at all 
sites. The time interval for linear splines was 60 minutes. The 
data were carefully analyzed for the presence of any remain- 
ing errors in group delay ambiguity resolution and for outliers. 
About 1 % of the data were deselected, primarily because the 
fringe fitting algorithm picked a secondary maximum of the 
delay resolution function. If only two observations of a source 
are used in a solution, errors in group delay ambiguity reso- 
lution cannot be detected. However, if more observations are 
used in a solution, errors in group delay ambiguity resolution 
can result in a significant misfit. Therefore, analysis of obser- 
vations of sources with numerous misfits allowed us to detect 
group delay ambiguity resolution errors with at least 3-6 ob- 
servations. Being conservative, we set a minimum number 
of 8 observations needed to rule out the possibility of group 
delay ambiguity resolution errors. After outlier elimination 
we determined additive baseline-dependent corrections to the 
a priori weights in such a manner that the ratios of the sums 
of weighted squares of post-fit residuals to their mathematical 
expectation were close to unity. 

In the final solution all available astrometric and geodetic 
data from August 1979 to August 2004, including the three 
VCS3 sessions, were used in a single least square solution. 
Estimated parameters included positions and velocities of sta- 
tions, coordinates of all sources. Earth orientation parame- 
ters, clock functions and troposphere zenith path delays. The 
only differences in treatment of the VCS3 sessions with re- 
spect to other experiments were that we did not estimate tro- 
posphere gradients and the time interval for clocks and tro- 
posphere zenith path delays were 60 minutes instead of 20 
minutes. Since a few sources were detected in only one band, 
we made three solutions: using 1) ionosphere free linear com- 
binations of X-band and S-band wide-band group delays, 2) 
only X-band group delays and 3) only S-band group delays. 

Since the equation for time delay is invariant with respect to 
the group of coordinate transformations, observations them- 
selves do not determine source positions, but only a family 
of positions. In order to pick a specific element of this fam- 
ily, we imposed boundary conditions by applying net-rotation 
constraints on the positions of 212 sources listed as "defined" 
in the ICRF catalogue such that the set of new positions of 
these 212 sources did not have a net-rotation with respect to 
the set of old position from the ICRF. More det ails about this 
procedure can be found in Be aslev et a 

In order to evaluate position errors we made a special so- 
lution using only the three VCS3 sessions in which we esti- 
mated coordinates of all sources and applied no-net-rotation 
constraints on 20 sources common in the VCS3 and ICRF 
defining lists. Since 77 sources used in the VCS3 campaign 
were intensively observed earlier under various geodetic and 
astrometric programs, we used the differences between very 
accurate positions derived from analysis of these observations 
and positions from the VCS3 sessions as a measure of the up- 
per limit of VCS3 catalogue errors. We represent estimates 
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TABLE 1 

Additive corrections to formal uncertainties applied for derivation of errors 

OF source position in VCS3 CATALOGUE 



Declination 




X/S 






X-band 






S-band 




zone 


#pts 






#pts 






#pts 


ha 




(deg) 




(mas) 


(mas) 




(mas) 


(mas) 




(mas) 


(mas) 


( +20.0, +90.0 ) 


33 


0.2 


0.4 


96 


0.6 


0.5 


94 


5.1 


4.7 


( 0.0, +20.0 ) 


16 


0.3 


0.2 


77 


0.5 


0.6 


11 


4.2 


7.2 


( -20.0, 0.0 ) 


11 


0.2 


0.3 


64 


0.8 


1.5 


64 


7.0 


16.8 


(-30.0, -20.0) 


11 


0.2 


0.7 


44 


1.2 


1.8 


44 


13.0 


27.6 


(^5.0, -30.0) 


13 


0.3 


0.4 


76 


2.4 


7.3 


75 


40.0 


40.0 



of errors in right ascension (ctq) and in declination {as) of the 
VCS3 catalogue in the form: 

a„ = V(l-5-02 + (£„('5)/co* 5)2 



where uj^ and crj stand for the formal uncertainties de- 
rived by the errors propagation law from fringe amplitude 
signal to noise ratios, Ea and Es are additive reweighting 
parameters. The multiplicative factor 1.5 in (|2}, first found 
by Rvan et al. ( 1993), accounts for underestimated system- 
atic errors in the computation of uncertainties of wide-band 
group delays and affects uncertainties of estimates of all pa- 
rameters from VLBI, including source coordinates. We split 
the sources into five declination zones, and for each group 
we found the values of Ea and Es which made ratios of the 
weighted sum of squares of residuals to their mathematical 
expectation close to unity. Using this techniques we com- 
puted additive corrections to uncertainties of the estimates 
of X-band only and S-band only solutions by analyzing the 
differences between dual band and single band solutions. 
Reweighting parameters are presented in Table [2 The his- 
togram of source position errors is presented in Figure [l] 

For imaging purposes the data were re-processed with the 
global fringe fitting procedure (task "FRING" ) with SNR 
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Position error, mas 

Fig. 1 . — Histogram of semi-major error ellipse of position errors. The last 
bin shows errors exceeding 9 mas. 



cutoff 4 at each individual channel. Global, station-based 
fringing allows detection of weaker sources than the baseline- 
based fringing used for astrometric analysis. (We did not use 
global fringing for astrometric analysis, since as we found 
earlier, it gives less reliable group delays than baseline-based 
fringing). A bandpass correction was then determined with 
the task "BRA^SS". After initial editing the fringe amplitude 
and phases were averaged across the channels. We used the 
AIPS task "CALIB" in order to perform the first step of the 
self-calibration procedure applying the point source model. 
The whole scan length was used to find the combined solu- 
tion for the channels with the SNR cutoff 3. Aft er that we 
imported data to the Caltech DIFMAP package ( Shepherd 
.1997,) . We have adapted the automated procedure origi- 
nally written by Greg Taylor (Pearson et al. 1994) for self- 
calibration and imaging of the VCS3 short snapshot observa- 
tions. 

For each detected source at each band, we computed two 
estimates of correlated flux density: one averaged over the 
entire image (total CLEAN flux density) and another aver- 
aged only over the spatial frequencies > 45MA at S-band 
and and > 170MA at X-band, which corresponds to the range 
[0.7, 1.0] of the maximum spatial frequencies of the VLBA. 
The latter estimate characterizes the flux density of unre- 
solved components. Errors of our estimates of total corre- 
lated flux density of sources brighter than 100 mJy are deter- 
mined mainly by the accur acy of amplitude caUbration, which 
for the VLBA, according to Wrobel and Ulvestad (2004), is at 
the level of 5% at 1-10 GHz. This estimate was confirmed by 
the comparison of the correlated flux density with the single- 
dish flux density measured with RATAN-600 in October 2003 
for several slowly varying sources without extended structure. 
The contribution of fringe amplitude errors is significant for 
sources with flux density < 100 mJy. The error of correlated 
flux density for this group of sources is about 10%. It should 
be noted that in the case of a source with asymmetric extended 
structure the estimate of the correlated flux density from un- 
resolved components may be biased due to a significant non- 
uniformity of uv coverage at long spacings. 

In total, 365 out of 450 sources were detected and yielded 
at least two good points for position determination. However, 
not all of these sources are suitable as phase referencing cal- 
ibrators or as targets for geodetic observations. We consider 
a source suitable as a calibrator if 1) the number of good X/S 
pairs of observations is 8 or greater in order to rule out the 
possibility of a group delay ambiguity resolution error; and 2) 
position error before reweighting is less than 5 mas. We found 
that usually 8 observations at the VLBA is sufficient to get 
positions with errors less than 5 mas. If the errors are signif- 
icantly greater, this means that the source was not detected at 
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TABLE 2 
The VCS3 catalogue 



Source name J2000.0 Coordinates Errors (mas) Con'elated flux density (Jy) 

8.6 GHz 2.3 GHz 



Class 


TVS 


lAU 


Right ascension 


Declination 


Aa 


A<5 


Con- 


#Obs 


Total 


Unres 


Total 


Unres 


Band 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


es) 


(9) 


(10) 


(11) 


(12) 


(13) 


(14) 


C 


0001-120 


J0004-1148 


00 04 04.914997 


-11 48 58.38564 


0.24 


0.40 


-0.069 


20 


0.62 




0.76 




X/S 




0003+340 


J0006+3422 


00 06 07.382435 


+34 22 20.41138 


1.83 


4.46 


-0.233 


5 


0.06 








X 


C 


0009+655 


JOO 12+6551 


00 12 37.671095 


+65 51 10.82363 


9.21 


4.00 


0.040 


10 


0.12 




0.69 




X/S 




0017+296 


J0019+2956 


00 19 37.791440 


+29 56 01.93471 


19.34 


7.63 


-0.015 


14 






0.08 




s 




0021-084 


J0024-0811 


00 24 00.672746 


-08 11 10.05397 


4.75 


8.88 


-0.124 


7 


0.08 


0.06 


0.09 


0.06 


X/S 


c 


0026-015 


J0029-0113 


00 29 00.986049 


-01 13 41.76047 


0.74 


0.95 


0.349 


36 


0.20 


0.06 


0.19 


0.10 


X/S 


c 


0027-024 


J0030-0211 


00 30 31.823755 


-02 11 56.13361 


0.30 


0.57 


-0.243 


68 


0.25 


0.22 


0.18 


0.16 


X/S 


c 


0029-147 


J003 1-1426 


00 31 56.411853 


-14 26 19.34670 


0.93 


1.60 


-0.607 


36 


0.17 


< 0.06 


0.23 


0.13 


X/S 


c 


0042+186 


J0044+1855 


00 44 42.227898 


+18 55 05.03460 


0.76 


1.69 


0.591 


36 


0.09 




0.14 


0.11 


X/S 




0042-373 


J0045-3705 


00 45 12.065955 


-37 05 48.47087 


3.05 


7.43 


-0.066 


20 


0.20 


0.19 






X 


c 


0043+246 


J0046+2456 


00 46 07.825737 


+24 56 32.52438 


0.26 


0.46 


0.051 


90 


0.44 


0.33 


0.18 


0.15 


X/S 


c 


0046+063 


J0048+0640 


00 48 58.723150 


+06 40 06.47544 


0.67 


1.20 


0.414 


28 


0.13 


0.06 


0.24 


0.06 


X/S 


c 


0054+161 


J0056+1625 


00 56 55.294324 


+16 25 13.34088 


0.33 


0.29 


0.019 


90 


0.42 


0.40 


0.27 


0.26 


X/S 


c 


0055+060 


J0058+0620 


00 58 33.804481 


+06 20 06.07358 


0.68 


1.31 


-0.221 


19 


0.17 


0.15 


0.09 


< 0.06 


X/S 



Note. — Tablel2lis presented in its entirety in the electronic edition of the Astronomical Journal. A portion is shown here for guidance regarding its form and contents. Units of 
right ascension are hours, minutes and seconds, units of declination are degrees, minutes and seconds. 



long baselines, and therefore, we do not possess information 
on whether the source can be detected when used as a phase 
referencing calibrator Only 23 1 sources satisfy this calibrator 
criteria. It should be mentioned that our criterion for suitabil- 
ity a source as a phase calibrator is rather conservative, and 
sources which fail this criterion may still be useful for some 
applications. 

Five detected sources turned out to be known gravitational 
lenses. Since the fringe amplitude may have several minima 
and maxima during integration time for such source, our tech- 
nique for fringe searching is not applicable. We excluded 
these gravitational lenses from the analysis, since they are cer- 
tainly not suitable as calibrators. 

4. THE VCS3 CATALOGUE 

The VCS3 catalogue is listed in Table|2] The first column 
gives source class: "C" if the source can be used as a calibra- 
tor, "-" if it cannot. The second and third columns give IVS 
source name (B1950 notation), and lAU name (J2000 nota- 
tion). The fourth and fifth columns give source coordinates at 
the J2000.0 epoch. Columns 161 and 111 give inflated source 
position uncertainties in right ascension and declination in 
mas, and column /8/ gives the correlation coefficient between 
the errors in right ascension and declination. The number of 
group delays used for position determination is listed in col- 
umn 191. The next two columns give the estimates of cor- 
related flux density in Jansky at X-band: total CLEAN flux 
density (column /lO/) and flux density of unresolved compo- 
nents (column /1 1/). Similar information for S-band is listed 
in columns /1 2/ and /1 3/. For some sources only estimates 
of the total VLBA flux density are presented, because either 
no data at long baseline were collected or the data at these 
baselines were unreliable. Column /14/ gives data type used 
for position estimation: X/S stands for ionosphere-free linear 
combination of X and S wide-band group delays; X stands 
for X-band only group delays; and S stands for S-band only 
group delays. Some sources which yielded less than 8 pairs of 
X and S band group delay observables had 2 or more observa- 
tions at X and/or S band observations. For these sources either 
X-band or S-band only estimates of coordinates are listed in 



the VCS3 catalogue, whichever uncertainty is less. 

In addition to this table, the html version 
of th is catalogue is posted on the Web at 
"http : / / gemini . gsfc.nasa. gov/vcs3" For each 
source it has 8 links: to a pair of postscript maps of the source 
at X and S-band; to a pair of plots of correlated flux density 
as a function of the length of the baseline projection to the 
source plane; to a pair of fits files with naturally weighted 
CLEAN images; and to a pair of fits files with calibrated 
uv data. The coordinates and the plots are also accessi- 
ble from the NRAO VLBA Calibration Search web-page 
http : / / magnolia .nrao.edu/ vlba_calib/ 

A sample of data used for candidate selection as well as 
some typical imaging results are shown in Figure|2l In row /a/ 
we present the broad-band spectra from RATAN-600 obser- 
vations in 2003 and from some old data taken from the NVSS, 
87GB (Gregory & Condon 1991) and the Radio Survey of the 
Galactic Plane at 1 1cm (Furst et al. 1990) catalogues. In rows 
Ihl and /d/ the dependence of visibility function amplitude 
versus projected baseline at S and X bands, respectively, is 
shown. Each point is a coherent average over the scan dura- 
tion and over the band. Naturally weighted CLEAN images 
at S and X bands are presented in rows Id and /e/"*. The 
sources in this figure are placed such that the contribution 
of the compact component in total emission incr eases from 
left to right. According to Kovale v et alJ (|2002'), we inter- 
pret the broad-band spectra of all selected sources as a sum 
of two constituents: emission from an optically thin extended 
component with steep spectrum which usually increases with 
decreasing frequency up to hundreds of MHz, and emission 
from a compact component with self-absorption which has a 
maximum at a frequency higher than 1 GHz. The spectrum 
of J0012+6551 has spectral index -0.5. No contribution to 
emission from the compact component is seen, and the corre- 
lated flux density even at the shortest baseline projections is 
several times less than the total flux density. See discussion 
of angular dim ensions of compact steep spectrum sources in 
lO'Dejfigga) for more details. The source is below the de- 

^ The scale of the J0012+6551 map is increased with respect to other maps. 
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Fig. 2. — a) Broad-band spectra from RATAN-600 observations in 2003 plus data from CATS, b) Dependence of the correlated flux density at S band 
(2.3 GHz) versus projected spacings. A dotted line shows the detection limit, c) Naturally weighted CLEAN images at S band. The lowest contour levels (2 
steps) on images are plotted at "clev" levels [Jy/beam], the peak brightness — "max" values [Jy/beam]. The dashed contours indicate negative flux. The beam is 
shown in the bottom left comer of the images, d) Dependence of the correlated flux density at X band (8.6 GHz) versus projected spacings. e) Naturally weighted 
CLEAN images at X band. 
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tection Hmit at both bands at spacings > 20 MA. At X-band 
the map shows only noise. The spectrum of J0909+4253 has 
the spectral index -1.1 at frequencies less than 2 GHz and 
it flattens at higher frequencies. This flattening is caused by 
the contribution of the compact component. When the con- 
tribution of the extended component is subtracted from the 
spectrum, the residual spectrum grows with increasing the 
frequency, and our measurements of the correlated flux den- 
sity confirm it. The remaining two sources have flat and in- 
verted spectra, which is an indication of dominance of the 
compact constituent in total emission. Again our VLBA ob- 
servations confirm it. Correlated flux densities at short spac- 
ings for J0056H-1625 and some other sources are even higher 
than the total flux density observed at RATAN half a year be- 
fore the VLBA measurements due to the strong variability of 
the object. Analysis of this sample illustrates the validity of 
our approach for candidate selection. 

5. CONCLUSION 

The VCS3 Survey has added 360 new sources not previ- 
ously observed with VLBl. Of these, 231 sources are suit- 
able as phase referencing calibrators and as target sources for 
geodetic applications. A careful strategy of source selection 
resulted in a very high detection rate; more than 80% of the 
sources in our candidates list were detected. 

This campaign contributed considerably in reducing the ar- 
eas with low density of known calibrators. If only sources 
from the ICRF-Ext2 catalogue are taken into account, 48.1% 



of the sky area above declination -45° has no calibrator 
within a 4° radius disk in any given direction. Combining 
the ICRF-Ext2 and the VCS 1 catalogues, this area is reduced 
to 13.7%. Adding the VCS2 catalogue reduces the area to 
10.5%. After the VCS3 campaign this area is only 5.3% of 
the sky above declination -45°. Figure |3l shows the area with 
insufficient calibrators density before and after the VCS3 ob- 
serving campaign. 
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